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Abstract 



The determination of the conductivity of a deterministic or stochastic 
classical system coupled to reservoirs at its ends can in general be mapped 
onto the problem of computing the stiffness (the 'energy' cost of twisting the 
boundaries) of a quantum-like system. The nature of the coupling to the 
reservoirs determines the details of the mechanical coupling of the torque at 
the ends. 

1 Introduction 

The transport properties of physical systems reserve many surprises, particularly in 
low dimensions, and despite many decades of efforts, a general theory is surprisingly 
not yet available [U [2] . This is true even at the classical level. 

The conductivity of low dimensional systems is often anomalous, with transport 
coefficients diverging with the system size. In such cases, the nature of the contact 
with the reservoirs at the boundaries becomes an issue. For example, the thermal 
conductivity of a finite chain becomes zero not only in the obvious case in which 
the contact is bad, but also in the limit when it is too good [21 E]- A conductivity 
computed on the basis of a closed chain without reservoirs can clearly not take 
into account these effects, and in any concrete physical realization of transport with 
anomalous properties, one has to consider the system as composed of both the bulk 
and the bath. 

Spin and charge currents in closed quantum chains have been related to the 
corresponding stiffnesses [HE]. This relation between a transport and an equilibrium 
property has been very fruitful. More recently, the question of a quantum open chain 
without baths has been addressed [6]. 
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In this paper we follow a different path. We shall establish a general correspon- 
dence between i) the transport properties of any classical (stochastic or determinis- 
tic) system in contact with reservoirs at its boundary, and ii) the stiffness (or helicity 
modulus [?]) of a zero-temperature system that is perturbed at the ends by a twist 
applied through two elastic 'handles' (see Figure 1). In this mapping, the details 
of the coupling between these handles and the system are important and reflect the 
nature of the coupling between reservoirs and the lattice in the original setting. 

As we shall see, the correspondence we discuss here works at a different level from 
the one of Kohn [3] and Shastry-Sutherland [3] . Here, the stiffness one calculates is 
that of the evolution operator [5] and not of the energy itself 0. 

The plan of the paper is the following. In the next section we introduce our set- 
ting, recall the definition of conductivity and of stiffness and explain our strategy. 
We will find convenient to use the operator bra-ket notation by which the probability 
distribution of Eq. (0Q) is associated to a Hilbert space quantum-like "state"; this is 
described in Section 3, together with two examples that will be repeatedly used in 
the paper. Armed with this preliminaries, we first show in Section 4 that the bulk 
current in a transport model subject to a gradient in the boundary conditions can 
be expressed in terms of a 'helical' operator. We then use linear response theory to 
rederive in Section 5 two equivalent expressions of the finite volume conductivity: 
the first coincides with the standard Green-Kubo formula (i.e., time integral of the 
bulk current autocorrelation function); the second one involves the autocorrelation 
function of an operator which depends only on the boundaries. This second ex- 
pression allows us to show that the conductivity is proportional to the stiffness of 
the evolution operator when a twist is exerted at the boundaries. This is shown in 
Section 6, which is then followed by conclusions. 



2 Transport model 

Consider an extended system made of iV components whose temporal evolution is 
either deterministic (Hamiltonian) or stochastic and whose boundaries are coupled 
to reservoirs. A convenient description of the system is given by the evolution 
(Fokker-Planck, Liouville,...) equation 

d 

— /j,(x,t) = -Hfj,(x,t) , 

fi(x,t) = e~ tH fi(x,t = 0) (1) 

where 

• fi(x,t) is the probability distribution in a point x of the phase/configurations 
space at time t. 

1 And in the quantum generalization of the present work, the twist will act on the closed time 
path action, involving twice as many fields, and not on the original Hamiltonian. 
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Figure 1: Stiffness of a bar twisted by elastic handles. 



• H is the evolution operator which is assumed of the form 

H = H B + \(H L (a L ) + H R (a R )) . (2) 

Here Hb denotes the bulk evolution, whereas Hl(ocl) and Hr^r), which depend on 
a parameter a, describe the interaction with reservoirs connected at the left/right 
boundaries, respectively. We have made explicit the strength of the coupling to the 
reservoirs with the parameter A. 




2.1 Conductivity 

We consider situations when the bulk evolution has a constant of motion E(x) - 
typically the energy or the mass. We shall assume E{x) is spatially local: it is the 
sum of N local contributions 

N 

E = J2 E *> (3) 

4 = 1 

each Ei a function of the site i or its near-neighbors. 

A probability distribution function that is concentrated on an energy (or any 
other conserved quantity) shell 5[E{x) — E ] satisfies: 

E(x)fi(x,t) = E of i(x,t) (4) 

so that the Hilbert space breaks into subspaces corresponding to different values of 
the conserved quantity. If the stochastic evolution Hb is conservative, there are no 
transitions between shells, so Hb is of block form, each one corresponding to a value 
of E Q . This in turn implies that E commutes with the bulk evolution operator: 

[E,H B ] = 0. (5) 

The reservoirs violate the conservation of E 

[E,H L }^0 ; [E,H R ]?Q, (6) 
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and, in addition, if the values of the parameter «£ ^ a R are different, a current 
is established, a stationary state eventually sets in and one has transport of the 
quantity E from one boundary to the other. Examples are: 

- when the parameter is the temperature, a = T, one has transport of energy; 

- when the parameter is a particle density, a = p, one has transport of mass. 

If an — oil is small, in the stationary regime the averaged total current across 
all links {J) is linear: 

{J) ~ K N (a R - a L ) (7) 
which defines the conductivity kn for a system of size N . 

2.2 Stiffness 

In order to recall the definition of the stiffness let us consider a simple example totally 
unrelated to the context of transport models. Suppose we are given a quantum 
Hamiltonian 

H = H B + X(H L + H R ) (8) 

with a bulk part 

N q2 N-1 

H B = - J2 h2 fT2 ~ J2 cos(x t+1 - Xi) (9) 
i=i ax i i=i 

which clearly has symmetry with respect to simultaneous shift of all angles. The 
expectation of the global angle is fixed by two boundary terms: 

H L = - cos(rri) ; H R = - cos(x N ) (10) 

which impose an optimal average profile (x^ = for alH G {1, . . . , N}. If we wish to 
calculate the stiffness, taking into account the specific coupling, we do the following: 
we compute first the lowest eigenvalue e° of Hb + \(Hl + Hr). Next, we twist the 
two 'handles' 

^twisted = _ + . H R w ^ d = - C0S (x N -9) (11) 

and calculate the new lowest eigenvalue e 9 . In the limit where the coupling A — > 
then, of course, imposing the twist has no effect. In the opposite limit A — > 00 
then x\ = —9 and xn = are strictly imposed. For a finite value of the coupling 
< A < 00 there is a competition between the bulk lattice which prefers a flat 
profile and the boundaries which force a gradient. In the framework of elasticity 
theory, the stiffness measures the cost of twisting and is defined from 

£< - £ °^(A^T)< 2e ' 2 - < 12 > 

valid to first order in 9 2 . The factor N — 1 assures a good scaling for large iV in a 
system with normal elasticity. 
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To define the stiffness of the transport model, which includes bulk and reservoirs, 
we proceed in a similar manner. We think of the evolution operator H in formula 
([T]) as a quantum-like, albeit non- Hermit ian, Hamiltonian. We can interpret (0) as 
an invariance of the bulk with respect to 'rotations' of angle 9 generated by a group 



e i9E 



The boundary couplings break this invariance (cfr. (jHJ)), we can use them as 
'handles' to impose a 'twist' by applying the transformation in opposite directions 
at the ends: 

^twisted = e i9E HRe -i6E . H twisted = e -i0E HLe i0E (13) 

The (zero-temperature) stiffness a of the handle+bulk system is defined as in (|12[) 
from the increase in the lowest eigenvalue of H tw[sted (6) = H^ isted +H B + fljwisted 
for small 9. 

The results we shall show in what follows is that the conductivity k is, up to a 
trivial factor, equal to the stiffness a. In general, anomalous conductivity amounts 
to the system becoming rigid, due to long-range correlations. When the coupling to 
the bath is weak, the conductivity is small: the stiffness is also weak because the 
handles turn without affecting the system. In the opposite case, when the coupling 
is too strong, the conductivity of a system with anomalous diffusion may go to zero. 
From the point of view of elasticity, what happens is that the handles twist strongly 
the first and last sites of the chain, and it is more favorable for a stiff system to 
concentrate all the twist between the first two (and between the last two) sites. 



3 Bracket Notation 

A convenient way to write the evolution equation ([I]) is provided by the bra-ket 
formalism. The probability distribution at time t is encoded in the state \ip(t)), 
namely 

\i/j(t)) = [ dxfi(x,t)\x) , (14) 



where \x) denotes a vector which together with its transposed (x\ form a complete 
basis of a Hilbert space and its dual, that is 

(x\x') = S(x - x') . (15) 

It immediately follows that 

(a#(t)> = MM) (16) 

and the evolution equation ([T]) takes the form of a Schrodinger equation with imag- 
inary time 

j t \m) = -H\m) ■ (17) 

To compute expectation we introduce the flat state 

dx (x\ (18) 
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which is such that 

(-\x) = 1 . (19) 

Then for any observable A we have that its expectation value at time t can be 
written as 

(A(t)) = J dxfi(x,t)A(x) = (-\A\if>(t)) . (20) 
Conservation of probability implies that 

(-\H = 0. (21) 

The stationary state \ip) satisfies 

H\ip) = . (22) 

For an isolated system, namely A = 0, the invariant measure are given by any 
arbitrary function of E, as it is immediately seen from Eq. (JSJ)- For example the 
microcanonical ensemble is given by the uniform measure 

dx5(E(x) -Ne)\x) . (23) 

In the presence of reservoirs, namely A ^ 0, the stationary state is in general not 
known. The boundaries operators Hl/r representing the action of the reservoirs are 
chosen such that, when the bath parameters are equal = ccr = a the unique 
invariant measure is given by the equilibrium Boltzman-Gibbs measure 

p e ~aE(x) 

|V)a = J dx^^\x) , (24) 

where Z is the normalizing partition function 

In the following we will denote the equilibrium Boltzman-Gibbs state \ip)a with 
| a). In case the system has a discrete configuration space, integrals in all the previous 
formulas are replaced by discrete sums and Dirac delta functions are replaced by 
Kronecker delta functions. 

The relation between a and a is easy to obtain for each problem using the fact 
that the equilibrium state \a), see Eq. 0211), is annihilated by the boundary 'bath' 
operators Hl/r(q) 

(H L/R (a))\a) = 0, (25) 

This is because, by assumption, the bath terms leave the corresponding equilibrium 
structure stationary. 



Example A: Hamiltonian systems 

A large class of systems is represented by a model composed of iV point-like 
particles interacting with their nearest neighbours. This is described by the Hamil- 
tonian 

N „2 N-l 

fi = E|+E%ri ( 26 ) 

8=1 1=1 
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The evolution equation ([!]) holds with x = (q%, . . . , qN,Pi, ■ ■ ■ ,Pn) a point in phase 
space. The bulk evolution operator is the Liouville operator 

N d dV d 

= = < 27 > 

where { , } denotes the Poisson brackets, that is for any functions / and g 

If \ = Y——-—— (28) 
^ % dpi dpi dqi ' 

The constant of motion for the bulk evolution is the total bulk energy E. Indeed 
Hb obviously commutes with the operator that multiplies by E: 

H B (Ef) = {Ef, E} = E{f, E} + f{E, E} = E{f, E} = E(H B f) Wf . (29) 

To model the interaction with the boundaries, stochastic heat reservoirs can be taken 
as Ornstein-Uhlenbeck processes whose variances T^m specify the temperature of 
the bath. This is represented by boundaries evolution operators 

- H ^ = wX T ^ + ^' < 30) 

- ffR(r * )= ^Nb +p ~) • (3i) 

When the two thermal reservoirs have the same temperature Tl = Tr = T, the 
stationary state \ip)p is the Boltzman-Gibbs equilibrium measure 

e -f3E(x) 

dx^—\x) , (32) 

where the inverse temperature (3 = 1/T is such that 

(E) = (-\E\P)=T . (33) 



Example B: Simple symmetric exclusion process 

Another class of models we consider are continuous time Markov processes. With 
the restriction of a finite configuration space S, the process is specified by assigning 
the rates of transition w(a', a) for jumping from a configuration o G S to a config- 
uration a' G S. The master equation for fj,(a,t), the probability distribution of a 
configuration a at time t, then reads 

= £ (w(a, a'Ma', t) - w(a', a)»(a, t)) = -H»(a, t) . (34) 
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To be more definite we will work with the simple symmetric exclusion process which 
can be expressed - in a operator formalism - as an SU(2) ferromagnet of spin 1/2 [§]. 
This corresponds to the stochastic process on the lattice {1, . . . , N} where particles 
jumps at rate 1 to one of their neighbours and each site can accomodate at most 1 
particle per site. Configurations n e {0, 1}^ are then identified with ket states 

|n) = | ni, . . .,njv) = ®f =1 \rii) , (35) 

which specify the occupation numbers of each sites, namely rij e {0, 1}. The bulk 
evolution is given by the transition rates 

w(n i+1 '\n) = -(n i ' <+1 \H B \n) = (l-n i )n i+1 

w{n Li+ \n) = -(rj> i+l \H B \n) = m(l - n i+1 ) . (36) 

where n lJ is the configuration which is obtained from the configuration n by remov- 
ing a particle in % and adding it in j. 

In operator notation, this is generated by [9] 

N-i / y 



i=i 

where the 5"s operators act as 



Hb — X) ( + &i + 2SiSi +1 - -J , (37) 



S?\ n i) = (1 -ni)\rii + 1) 
S~\rii) = rijlrii-1) 

Si\ni) = (rij - - ) \rii) . (3^ 



and satisfy the SU(2) algebra 



[slsf] = ±st 

[Sr,S+] = -2fi?. (39) 



In this stochastic system the constant of motion - in the absence of reservoirs - 
is the total number of particles 

E = J2 (5 + ~) ■ (40) 

This is obviously a conserved quantity for a particle jump process, 

[E,H B }=0, (41) 



as can be immediately checked using the commutation relations (1391) . 

When the system is coupled to reservoirs at different chemical potential, the 
boundaries can inject or absorb particle from the system. We assume that particles 
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are injected at site 1 at rate p L and they are removed from site 1 at rate (1 — p L ); 
in the same way particles are injected at site N at rate p R ) and they are removed 
from site N at rate (1 — p R ). The boundaries evolution operators then read 



- H L ( PL ) = PL (s+ + Sj - I) + (1 - p L ) (S; - S[ - ±) (42) 

- H R { PR ) = PR (s+ + S° N -l) + (l- p R ) (s N -S° N -l) (43) 



The constants pi and p R are to interpreted as the two densities the reservoirs would 
impose if acting separately. When the two reservoirs have the same densities pi = 
Pr = P, the stationary state \ip) u is the Boltzman-Gibbs equilibrium measure 



vE 

H = EV |n) ' (44) 

{n} 



where the chemical potential v = In (j^j is such that 

(E) = (-\E\u)=p. (45) 

As it is immediately checked, in this case the equilibrium Gibbs measure is a 
Bernoulli product measure with parameter p = e v /(l + e u ) 

\P) = ®li\Ph = |1>< + (1 - P) |0><) • (46) 



4 Bulk current and 'helical' operator 

We wish to study the transport of the quantity E in the presence of reservoirs. Our 
first step is to express the current in the bulk in terms of a 'helical' operator. 

The current which passes through the system is defined via the continuity equa- 
tion 

f + VJ = 0. (47) 

Since E is the sum of iV local contributions, see Eq.(j2J), we can write a local conti- 
nuity equation as 

dE 

-^ = -(Js,s+l-Js-l,s), s = l,...,N. (48) 

where J s , s +i is the current from site s to site s + 1, and the total bulk current is 
defined as 

JV-l 

J = E Js,s+i ■ (49) 

8=1 
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For a stochastic system, J s>s+ i may depend explicitely on the random noise. We 
shall consider a set of operators J s>s +i such that they coincide with the average of 
Js, s +i over random realisations 

(Js, s+ i(t)) noisc = (-1 J.,,+1 W)> ; (J(t)) noisc = (-\J\y(t)) (50) 

Defining the helical operator 

N 

A = Y,sE s . (51) 

s=l 

J can be expressed as the commutator between it and the bulk evolution operator, 
that is 

J=[A,H B \ (52) 
Indeed, for an open chain (namely A = 0) we have 

1*0)) = • (53) 



Making use of the continuity equation (|48p . we then have 

d A N dF N 

(-\[A,H B ]Mt)) = H^W)> = £«H^W)> = 5>U-i, s - ■ 



(54) 

Since J7b,i = Jn,n+i = for an open chain: 

(-|[A,^]|M/(t)) = f: S (-|J s _ M |vI/(t))-^ 1 s (-|J SiS+1 |vl/( t )), (55) 

s=2 s=l 

and making the change s — > s + 1 in the first summation we find 

N-l 

(-\[A,H B ]Mt)) = = W = HW)> • (56) 

8=1 

Remark 1: The correspondence between J and J is via expectation values. 
Note however that expectations of products yield extra terms in stochastic systems: 
(J 2 ) = (— | J 2 |^/) + extra trems. This is the origin of extra terms in the formulas 
below. 

Remark: 2 In the definition of the helical operator one can always add a term 
whose total derivative with respect to time is zero. Later, to make the role of the 
boundaries more symmetric, we will consider 

A N + 1 

A = J2sE s - 7j~E (57) 
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Example A: Hamiltonian system 

For the system described by the Hamiltonian (T2"6"]) the current from site s to site 
s + 1 is given by (see [2]) 

Js,t+i = ~\(Ps + Ps+i)V\q s+1 - q s ) (58) 
Defining the local operators E s as 

E * = 4 + l(v{qs + i-q s ) + V{q s -q s _ 1 )) s = 2, . . . , N - 1 



E N = V ^ + \v{q N -q N ^) (59) 
one can check that (1521) holds with 



8=1 Z 

N v 2 N ~ 1 N + 1 

= £^ + £(2s + i)n? s+ i-^)-^£ (60) 

g=l A S =l / 



Example B: Simple symmetric exclusion process 

For the system described by the evolution operator (1371) the current from site s 
to site s + 1 is 

<^s,s+i = S s S s+1 — S s S s+1 (61) 
Equation (1321) now holds with 

4.1 A useful general identity 

The boundary terms Hl/r{ol) break the conservation law for the quantity E. They 
depend on a parameter a (e.g. the temperature in the Hamiltonian system example, 
the chemical potential in the symmetric exclusion process example) and impose a 
unique stationary equilibrium Boltzmann Gibbs state \ip)a = \oc)- We wish to find 
a convenient way to express the response of the boundary operators H L m(a) when 
their parameter a is varied. We claim the following identities hold: 

= c Q [H L/R (a),E] \a) , (63) 
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where the constant c a is given by 



da 
da 



(64) 



This can be proved as follows. We use the fact that the Gibbs-Boltzmann distri- 
bution is annihilated by the boundary operators (cfr Eq. (1251) ). and that, because 
of the Gibbs-Boltzmann form, it satisfies: 
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-(E-(E))\a) 



(65) 



where (E) = (—\E\a). Computing the derivative of Eq. (l25p with respect to a, this 
in turn gives: 







d_ 

da 

da 



\ d 
—H L/R (a) I | a) + H L/R (a) — \a) 

z-H L/R (a)\ | a) - H L/R {a){E - (E)) \a) 
d \ 

—H L/R (a)\ | a) - [H L/R (a),E\ \a) 



(66) 



where we have used 



Eq. 



then follows. 



Example A: Hamiltonian system 

Eq. ( 1631) holds with a = T 1 a = f3=k and with ct 



■ T 2- 



Example B: Simple symmetric exclusion process 

Eq. fl6"3j) holds with a = p, a = —v = In (j-^j and with c p 



i 



p(i-p) ' 



5 Linear response theory 

The conductivity is obtained by calculating the average current in the presence of a 
small mismatch in the parameters of the reservoirs evolution operators. If the left 
reservoir is working with a parameter a^ = a — 5a and the right reservoir is working 
with a parameter a R = a + 5a, the evolution equation 

j t \m) = -mm) ^) 

can be solved in a linear response regime. Expanding to first order in 5a we find 

H = H + 5a X AH + o(5a) (68) 

with 

H = H B + \{H R {a) + H L (a)) (69) 
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AH= {l H *-i H i) (70) 

and 

\iji(t)) = \a) + 6a \Aijj(t)) + o{6a) (71) 
with | a) is the stationary state for the unperturbed problem, that is 

j t \a) = -H \&) = (72) 

The solution is 

\Aip(t)) = -X f dt' e- {t - t,)H AH\a) (73) 
Jo 

To leading order the average value of the total bulk current will be 

(j(t)) = (-\j\m) 

= (-\J\a) - 6a\ [ dt' (-\Je-^ H AH\a) + o{6a) (74) 
Jo 

The conductivity for a system of size N is defined as 

- lim ™ (75) 
From Eq. (174p the following formula for the conductivity is then deduced: 

k n = lim -- C dt' (-Ue^-^AHla) (76) 

t^co 2 JO 

Thanks to the previous established (1631) . we have: 

A#|a) = -ca H'\a) (77) 

where we have defined 

H'=[H L -H R ,E}. (78) 
Substituting in ( 1761) . we then have 

k n = lim ^ /V l-\Je-^ H H'\a) (79) 
5.1 Green-Kubo formula in the bulk 

The Green-Kubo formula for the conductivity usually involves the current tempo- 
ral autocorrelation function [TU] . To see that formula ( 1791) is indeed the standard 
Green-Kubo formula we proceed further by using some properties of the reservoirs. 
Specifically, we use the assumption that the reservoir on the left (resp. on the right) 
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depends only from the phase/configuration variables of the first (resp. last) site. 
This implies that 



N 



[A,H L + H R ] 



1 

N 



■N N+1 N 

J2sE s I-J2E S ,H L + H R 

.8=1 1 8=1 



1 




N 


[(" 


N - 


- 1 


2N 


N - 


- 1 


2N 



'I -N 



1 

+ iV 



vV-1 



[H l -H Rj E] 
H' 



JO) 



On the other hand, we also have 



i[Afl-Ji + ffL] = ^[Aff--fffl] 



Putting (IHUj) and (IHTj) together we find 



A(iV 



52) 



The standard Green-Kubo formula is obtained by substituting H' on the right of 
SJ) by the expression found in Eq. fl82|) . We have 



lim 

t-»oo jV — 1 



dt' (-|Je- ( '-*' )H J|«) 



-lim T7~~r r^(-|J e - ( '-'' )H [i/,A]|a) (83) 

t^oo iV — 1 JO 

The term involving the commutator [H, A] can be further simplified as follows 
lim C dt' {-\Je- {t - t,)H HA\a) = lim /' dt' A.(-\J e -^ H A\a) 



(-\JA\a) 



54) 



where we have used that H\a) = and we have assumed that the J and A have 
vanishing correlation for very large times. Moreover, by using ( 1321) and the fact that 
Hb\o) = { — \Hb = 0, one has that 



(-\JA\a) = (-\[A,H B ]A\a) 
= (~\AH B A\&) 

= l(-\[[AH B },A}\a) 



55) 



Finally we get 



Km = — lim — 

JV t—oo j\T 



^- fdt {-\Je-W H J\a) 
— 1 Jo 



2(N-r 



(-\[[A,H B ],A]\a) (86) 
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The second term comes from the fact, already mentioned above, that ( — | Je~(* - * ^ J\ct) 
is, in stochastic systems, the current-current correlation function only up to an equal- 
time extra term. We distinguish two cases: 

• A = 

In this case H\ = H B and then the current autocorrelation function can be 
treated as we just did for the extra term: 

lim f dt' {-\Je-^ HB J\a) = lim /* dt' (-1 Je- {t - t>)HB H B A\a) 

= lim tdt' -^(-\Je-^ HB A\a) 

t^ooJo dt'^ ' 1 ' 

= (-\JA\a) 

= ~(-\[[A,H B ],A]\a). (87) 

This implies that the two terms in Eq. ( 1861) cancel each others and, consistently 
with the fact that there is no coupling to the reservoirs, we find = (see 
Ref. [6] for a discussion) 

• A ^ 

This time is different from zero and is, in general, the sum of two competing 
terms: the time integral of the current autocorrelation function and an extra- 
term which is purely due to randomness. 



Example A: Hamiltonian system 

The commutator in Eq. (1861) vanishes in the Hamiltonian case. Indeed we 
have: 

[A,H B ]f = A{E,f}-{E,Af} V/ (88) 

from which it follows that 

[[A,H B ],A]f = [A, H B ]Af — A[A, H B }f 

= A{E, Af} - {E, A 2 f} - A 2 {E, /} + A{E, Af} 

= (89) 

Inserting a = T and ct = — ^ into ( 1861) . the Green-Kubo formula for the 
conductivity then reads 

K - = fe 55(173!) f^t'(-\Je-"- l ' >H Jm (90) 



Example B: Simple symmetric exclusion process 
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In this case the commutator in Eq. (186]) is non-zero. We have 



N-l 



[A, H B } = J=J2 [S7S: +1 - StS; +1 ) (91) 



8=1 

and 

N-l 

p, h b },a} = -J2 (s;s+ +1 + s+s; +1 ) (92) 

8=1 

Recalling that for the SSEP with reservoirs having the same chemical potential 
the equilibrium state is given by a Bernoulli product measure, see Eq. (l4"6"j) . 
an immediate computation gives 

(-\[[A,H B ],A]\a) = -2p(l - p)(N - 1) (93) 

This implies the following expression for the conductivity: 

«» = -! + fe (iV-lMl-p) t d< HJe-^JM (94) 

As a final remark of this section, let us check that in the thermodynamical 
limit the correct value of the conductivity is recovered. In order to evaluate the 
contribution due to the current autocorrelation function we observe that in the 
thermodynamic limit N — > oo, whatever the value of A < oo, the boundaries 
will be negligible. This lead us to evaluate this term for the infinite volume 
system: 

lim lim — -L- r C dt! (-|J e -(*-*')tf J\ p ) = 

N^oot^oo (AT _ l)p(l _ p) J Q N 1 IF/ 

/ dt' < Jo,i(0)Jo,i(0 > 
J 



P(l - P) 



(95) 



where < ■ > denotes expectation with respect to the equilibrium state. We 
can then use the duality property for the model and specifically the following 
results: 

< (n x (t) - p) (n (0) - p) > = p(l- p) Pt (x) (96) 

where pt(x) is the probability that a continuous time simple symmetric random 
walk jumping left or right at rate 1, started at the origin at time zero, is found 
at site x at time t. Using duality we have 

< Jo,i(0) J ,i(t') > = < (m(0) - no(0)) M*') - 710(0) > 
= p(l-p)(-p t ,(-l) + 2p t ,(0)-p t (l)) 

= - p (l- p )± Pt ,(0) (97) 
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Putting together (1951) and (1971) we arrive to 

fete (iv^lWw) i' dt '<-l Je """'' )Hj l"> = -2(P«(0)-P»(0)) 

= 2 (98) 
Inserting this result in Eq. f[9~4l we finally find (as it should be!): 

k= lim k n = -1 + 2 = 1 (99) 

5.2 Green-Kubo formula for the boundaries 

Here we follow the opposite strategy of the previous section. Recall expression (ff9l) 
for the conductivity: 

K N = lim ^ t dt' (-\Je-^- t,)H H'\&) (100) 

t^oo 2 Jo 

This time we express J in terms of H 1 by inverting relation f)82p . that is: 

J = - KN ~ 1) H'-[H,A] (101) 
This yields the following expression 

KN = -lim ^(^-i) f'dt' (-\H'e-^ H H'\&) 

t^oo 4 JO 

_l im A^ /"V (-ir^Ale-^^'la) (102) 

t^oo 2 JO 

As in the previous case the term involving the commutator [H, A] can be simplified 
by using the fact that it is the integral of a time-derivative, and that the correlations 
of A and H' vanish at widely separated times: 

KN = -lim A2 ^~ 1 ) f d f {-\H'e-^ H H'\&) 

t^co 4 JO 

-^(-\AH'\a) (103) 

The extra term can be rearranged as follows. Recalling the definition of H', Eq. 
(I7H1) . and using the fact that Hi and H R annihilate the equilibrium measure we 
have 

(-\AH'\a) = (-\A[H L -H R ,E]\a) 

= (-\A(H L -H R )E\a) (104) 
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Now we use that H L (resp. H R ) commutes with all the E i with i ^ 1 (resp. i ^ N). 
This yields: 



(-\AH'\&) = -^-±{-\E(H L + H R )E]\a) 



2 

N- 1 



[£,tf L + if fl ]]|a) 



(105) 



and the final expression for the Green-Kubo formula for the boundaries read: 

X 2 c a (N-l) ft. 



K N 



— lim 

t—>oo 



[ alt' (-\H'e-^ H H'\a) 
Jo 



8 



(-\[E,[E,H L + H R ]]\a) 



(106) 



We will show below that the extra term is non-zero for both the Hamiltonian chain 
and the SEP and has actually the same value. 

Example A: Hamiltonian system 
The extra term is evaluated as 

(-\[E, [E,H L + H R }]\P) 



We have 



H£i#l£i|/3> 



(-\[E U [E 1: H L }}\/3) + (-\[E N , [E N ,H R ] 
^(-^HlE^P) - 2{-\[E N H R E N \P) (107) 



t(-\pI\P) 



rj-i2 



(108) 



where we have used the fact that annihilates the flat measure to the left and 

op l 

(T-t^ + pi) annihilates the equilibrium measure to the right, together with 



p\ d 
~2'dp[ 



-Pi 



(109) 



d 



, pi 
T dp~, +Pl ^ 



Analogously we find that 

(-\E N H R E N \(3) = T 2 
so that the final result for the extra term in the thermodynamic limit is: 



(110) 



(111) 



(112) 



18 



Example B: Simple symmetric exclusion process 
With 

Hl = -P (V + S° - ±) - (1 - p) (V - S? - (113) 



we have 

(-|[E,[E,^ L ]]|p) = (-![£?!, [^ffjiip) 

= (-\-pSt-{l-p)S^\p) (114) 

Because (5"^ + — |) and (S^ — S® — |) separately annihilate the flat measure, we 
can substitute by S^'s to get 

H[£,[£,tf L ]]|p> = - + (1 - ^) - i) b> 

= -2p(l-p) (115) 

The extra term is then: 

^4^ H [^,m*> = ^^ (116) 

6 Stiffness 

Let us now compute the stiffness. We have a 'quantum- like' Hamiltonian H = 
Hb+\{Hl+Hr) where the bulk term is "symmetric" with respect to transformations 
generated by E, since [Hb,E] = 0. Hl and H R are the boundary 'handles' that 
break the symmetry generated by E. In order to calculate the stiffness, we twist 
the boundaries, i.e. we apply the transformation in opposite directions: 

Hl — > H e L = e t9E H L e- WE H R ^H e R = e~ ieE H R e WE (117) 

As a consequence of the twist the spectrum of the transformed Hamiltonian H e = 
Hb + A (iff + H e R ) will be different from the one of the original Hamiltonian H . We 
define the stiffness a in terms of the lowest eigenvalue of H 6 : 

= ^ >n(He-»^( f ))) 

t^oo t V ' 

We define 

H d = H[ + H R (119) 
Developing up to order 9 2 we have 

AH = H e - H = e t6E H L e~ t6E + e- idE H R e idE - (H L - H R ) 

= -iQH' - ( L([E,[E,H L + H R }}) + ... 
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Application of time-dependent perturbation theory gives 
]n(<-|e-**|V>(t))) = + X J dl/(-\AH\&) (120) 



+ fdt' f dt'i-lAHe-V'-^AHla) + ..) 
2 Jo Jo J 



In (l-6 2 ^t(-\([E, [E,H L + H R ]])\a) 



p*. f dt' f dt"{-\H'e-^"-^ HB H'\a) + ..\ 
2 Jo Jo J 



We find for the increase in the lowest eigenvalue 

e ~ e ~ 2(iv3T) (2 ^ 

ln(<-|e-^(t)> 
= — hm — 

= -\im^{-\([E,[E,H L + H R ]})\a) 
- lim 9 2 X 2 C dt'(-\H'e~^ HB H'\a) 

t-*oc Jo 

Recalling the expression fll06p for the conductivity 

2KN - = -lim^ fdt' (-\H'e-^ H H'\a) 
c a (N-l) 2 io ^ 17 

- lim ^ (-![£, [22, 22 L + H R ]]\&) (121) 
one obtains the following relation between the conductivity and the stiffness: 

KN = C a (TN ■ (122) 

Stiffness and conductivity are proportional up to a trivial factor. 

7 Conclusions 

When a system has a conserved bulk quantity, its bulk evolution operator has a 
symmetry. Current transport induced by boundary 'reservoir' terms corresponds in 
all generality to a twist exerted applying the symmetry transformation in opposite 
senses to the boundary terms. The conductivity of the system is in this view the 
stiffness, or helicity modulus, associated with this operation. 

This mechanical analogy of transport can be taken further. For example, one 
easily understands that in an elastic system any local perturbation that couples with 
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torsion will have long-range effects: the analogy discussed here means that the same 
can be said of perturbations of systems with conserved quantities [TT] . 

In fact, one recognizes methods where auxiliary thermal baths are used, with 
their temperature fixed so that they exchange no current on average [12] as the 
usual symmetry-breaking fields, adjusted so that they exert no average force, of 
statistical mechanics. 

The derivation in this paper was made for classical and stochastic systems, in 
contact with reservoirs at the ends. The generalization to a quantum system with 
baths can be made through the Feyman- Vernon [13] formalism. It is not at this 
point clear to us how this may be related to the approach of Kohn [4] and Shastry 
and Sutherland [5], where the twist is applied to the Hamiltonian (rather than the 
evolution operator) of a closed chain. 
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